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Abstract 

We study properties of the energy minima obtained by quenching equilibrium configura- 
tions of the Sherrington-Kirkpatrick (SK) mean field spin glass. We measure the prob- 
ability distribution of the overlap among quenched configurations and the quenched 
energy, looking at the dependence on the starting equilibrium temperature, and per- 
forming a systematic analysis of finite size effects. 
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1 Introduction. 



The local energy minima properties are analyzed in several works on glass-forming liquids and 
they are allowing a better understanding of the behavior of these systems. The importance 
of potential energy landscape in the physics of super-cooled liquids was already pointed out 
by Goldstein |l[|. More recently Stillinger and Weber || formalized the idea that the multi- 
dimensional energy surface can be partitioned in a large number of local minima, so called 
Inherent Structures (IS), each one surrounded by its attraction basin. It is now clear || that 
the low-temperature dynamics (i.e. for temperatures below the Mode Coupling critical tem- 
perature Tmct H) can be subdivided into intra-basin motion and crossing of energy barriers 
by activated processes, taking place on a significantly longer time-scale. 

The system at equilibrium below Tmct is 'almost always' trapped in one of the basins ac- 
cessible at this temperature. The huge number Af oc exp(iVE) of these 'valleys', exponentially 
diverging with the system size N, suggested the scenario of an underlying thermodynamic 
transition due to an 'entropy crisis' at the Kauzmann temperature T% < Tmct where the 
configurational entropy S goes to zero ||, which was supported by recent analytical work ||- 
|| . By looking at IS one can evaluate S numerically ||-[l(J. Moreover, the IS energy turns 



out to be an interesting quantity for studying both the static and the dynamical behavior |5J- 
H| . Differences between fragile and strong glasses were recently proposed to be explainable 



within an energy landscape description, too ]TJ| . 

The outlined picture of glass-forming liquids is reminiscent of that characterizing general- 
ized mean field spin glass models like those involving p-spin interactions fl5j , which display a 
dynamical ergodicity breaking at the temperature Tp = Tmct (in this case the barriers 
between basins are infinite in the thermodynamic limit also for Tk < T < Tmct because of 
the mean field approximation), and a thermodynamic entropy driven transition at a lower 
temperature Tr, corresponding to a one step replica symmetry breaking (1RSB) scenario. 
For T < Tk one finds a non trivial probability distribution of the overlap between states 
P(q) = m5(q) + (1 — m)5(q — Qea) ■ Here qEA is the self-overlap of a state with itself, whereas 
different states are orthogonal and have mutual overlap zero. 

Several years ago, Kirkpatrick, Thirumalai and Wolynes [17] suggested that 1RSB spin 
glass models could be a paradigm of vitreous systems. The numerical study of out-of- 
equilibrium dynamics in glass-forming liquids gives intriguing results The mea- 

surement of P(q) among 'glassy states' is a subtle task [^, ^3], |24|, since one faces both 



the problem of thermalizing the system down to very low temperatures and that of avoiding 
possible crystalline minima whose basin of attraction could be non-negligible for small sys- 
tems. As recently proposed by Bhattacharya, Broderix, Kree and Zippelius |[23|| , to look at 
the Inherent Structures is helpful also from this point of view, since it allows a more precise 
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definition of the overlap and make easier to distinguish between glassy minima and crystalline 
or quasi- crystalline configurations. 

On the other hand, it is not completely clear a priori which kind of behavior one should 
expect for the P que n{Q) measured among energy minima obtained by quenching equilibrium 
configurations instead of among equilibrium configurations themselves. To our knowledge, 
such a quantity was never previously studied in a spin glass model (which is not surprising, 
since in this case equilibration is still feasible without difficulties for moderate system sizes). 
More generally, little is known about properties of 'Inherent Structures' in spin glasses, which 
(apart from analogies with glass-forming liquids) have their own interest even in the well 
understood Sherrington Kirkpatrick (SK) mean field model. In a previous work one 
of the authors studied the energy minima of the SK model obtained starting from random 
initial configurations (i.e. by quenching from infinite initial temperature) whereas recently 
Crisanti and Ritort |26|, [27] have performed both for a 1RSB spin glass and for the SK model 
a numerical analysis of the static and dynamical properties within the energy landscape 
description similar to the one proposed in ||, [13] for a glass-forming liquid, by looking in 
particular at the IS energy and at the configurational entropy. Their results confirm the close 
similarities between 1RSB spin glass and structural glass energy landscapes and the usefulness 
of this kind of approach, further suggesting a systematic analysis. 

In the present work, we perform an extensive numerical study of energy minima properties 
in the SK model, considering initial equilibrium configurations both in the high-temperature 
(paramagnetic) region and deep in the glassy phase. We look first of all at the behavior of the 
appropriately defined P que n(<l) which is compared with the corresponding (usual) equilibrium 
one. Then we extend the analysis to the overlap of the IS with the configurations from which 
they are obtained, which measure the (quite strong) correlations between them, and we study 
systematically finite size effects on the behavior of the IS energy. 



Model, Observables and Simulations. 



The Sherrington Kirkpatrick spin glass model p8] is described by the Hamiltonian 

N 

'Hj = Jij<Ji&j , (1) 
i<j=i 

where a, = ±1 are Ising spins, the sum runs over all pairs of spins and J^- are random 
independent variables with mean value J,.,- = and variance 1/N. We take = ±N~ 1 ^ 2 . 

This model is exactly solvable (since interactions have infinite range) and has a glassy 
phase with full replica symmetry breaking (FRSB). In case of zero magnetic field (the one 
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we consider here), taking into account the symmetry under inversion of the spins, the P(q) 
changes at the critical temperature Tc = 1 from a 5-function at q = (characteristic of 
the paramagnetic phase) to the FRSB two 5-functions in q = ±qEA with a non-zero plateau 
joining them. The transition is continuous also in the order parameter (at variance with 1RSB 
models), i.e. lim T ^ T - g£^(T) = 0, and there is no distinction between the dynamical and the 
static transition PSl, i.e. Tn = T, 



d — j-c- 

The SK model is particularly suitable for the kind of study we are interested in, since 
its behavior is well understood and, on the other hand, by using optimized Monte Carlo 
methods |3(J one is able to thermalize large system sizes down to low temperatures, which 
allows to study finite size effects systematically. We simulated N = 64, 128, 256, 512 and 
1024, averaging over 2048, 1024, 512, 384 and 192 different disorder realizations respectively. 
The program was multi-spin coded on different sites of the system (we store 64 spins in 
the same word) and we used Parallel Tempering (PT) |30], running simultaneously two 
independent sets of copies (replica) of the system for each sample. Up to 50 (for the two 
largest sizes) different temperatures between T min = 0.65 and T max = 3 were used, and we 
performed from 100.000 PT steps for the smallest value of N = 64 to 300.000 for the largest 
value of N = 1024. 

The PT acceptance for the exchange of nearest neighbor temperatures was never smaller 
than 0.6. In the second half of the run we computed the specific heat both as the deriva- 
tive of the energy density with respect to temperature c = d(e)/dT, and from fluctuations 
c = N((e 2 ) — (e) 2 )/T 2 , looking for compatibility of results. This means to compare one-time 
and two-time quantities respectively, which is an effective way for checking thermalization par- 
ticularly when using PT (note that in this case fluctuations involve different replicas evolving 
at the same temperature at different times during the run). Nevertheless we also divided the 
second half of the run in (four) equal intervals, checking that there were no evident differences 
in the values of the considered observables, P(q) and P que n{q) in particular. A further confirm 
for the system being well thermalized comes from the perfectly symmetric with respect to the 
exchange q — > —q probability distributions that we obtained. 

For each disorder realization and for each temperature, in the second half of the run, we 
saved 1500+1500 equilibrium configurations from the two independent sets of replica which 
were subsequently quenched by a zero-temperature dynamics. The observables were computed 
from these configurations and the corresponding energy minima, errors being evaluated from 
sample-to-sample fluctuations. One should note that both P(q) and P que n{q) are strongly not 
self-averaging in the glassy phase, wherefore it was necessary to average over a large number of 
samples even for large system sizes. The whole simulations would have taken about two years 
of CPU on a usual alpha-station, i.e. a few days when using 128 processors simultaneously 
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on Cray T3E (the code is easily parallelized with efficiency close to 1 by running a different 
disorder realization on each processor). 

A subtle point concerns the quenching procedure. We are considering even values for N, 
which means that the local field acting on each spin because of the other ones can never be 



zero. Moreover, it is known f25fl , from the analysis of properties of quenched configurations 
obtained starting from infinite temperature (i.e. random initial configuration), that different 
zero temperature dynamics give qualitatively identical results. The 'Greedy algorithm', where 
the spin corresponding to the largest energy decreasing is flipped at each step, seems to stop 
more frequently in local high energy minima than the 'Reluctant algorithm', where one flips 
the spin which gives the smallest energy decreasing (i.e., the contrary of the previous case). We 
choose to use an 'intermediate' zero-temperature dynamics, which is easy to implement too. 
At each step, a randomly chosen spin is suggested to flip and at least 20N steps are performed 
after the last successful one before stopping. The probability that the final configuration is 
not a local energy minimum (under single spin flip) is therefore ~ e~ 20 , practically negligible. 
As a last remark, it should be stressed that such a quenching procedure could possibly give 
energy minima which are not the 'nearest' IS to the starting equilibrium configurations, i.e. 
which are less correlated but not more correlated than in the analog glass-forming liquid 
case, where one generates ISs by following the path of steepest descent. 

Labeling by {t^} the spins belonging to two configurations, the overlap is defined as 

1 N 

G = -tf£<W- ( 2 ) 

i=l 

The spin glass order parameter, i.e. the equilibrium probability distribution of overlap among 
states P(q, T), is usually evaluated numerically as the histogram of the instantaneous overlap 
Q between two replica (with the same disorder configurations) which evolve simultaneously 
and independently at temperature T 



P(q,T) = Pj(q) = (5(q-Q)} , (3) 

where the thermal average (•) corresponds to average over time in the simulation whereas 
(•) stands for the average over realizations. Pj(q,T) can be equivalently measured from 
two given sets of equilibrium configurations belonging to the two replica, by considering the 
overlap of each configuration of one set with all the configurations of the other. In this work 
we evaluated Pj(q,T) both during the simulation and from the saved configurations, obtain- 
ing perfectly compatible results, which confirm that these configurations sample accurately 
enough the phase space. 
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We define the quenched probability distribution of the overlap as 



i Mis 

P quen (q,T) = £ 5{q-Q I8 ) , (4) 

•"is ia ,i b =i 

where the sum runs over the jV/s — 1500 energy minima obtained starting from the equi- 
librium configurations at temperature T for each of the two replica sets. This definition, 
analogous to the one introduced in for a Lennard- Jones glass-forming binary mixture, 



implies that we are weighting each IS with the Boltzmann factor of the corresponding basin 
at temperature T, as usual in numerical studies on super-cooled liquids. 



3 Results. 

3.1 The behavior of P qU en{q)- 

We present in [Fig. [l[ data for the equilibrium overlap distribution P(q) (on the left) and for 
the quenched overlap distribution P quen ((l) (on the right) at a temperature very close to the 
critical temperature T = 1.05 = 1.05 Tc, but still in the paramagnetic phase, for different 
system sizes. It was shown in |2H| that the P quen {q) obtained from infinite temperature 
configurations becomes more and more concentrated in q = for increasing N and it goes to 
a 5-function in the thermodynamic limit. We find the same qualitative behavior in the whole 
high-temperature region T > Tc = 1. It is crystal clear from [Fig. |IJ that there is no evidence 
for replica symmetry breaking in the probability distribution of overlap between IS reachable 
from equilibrium configurations at T ~ Tc and weighted with the Boltzmann factor of the 
corresponding basin at this temperature. 

In other words, we know that for T < Tc there is RSB, but we cannot detect it by 
looking at the probability distribution of the overlap obtained with a fast quench starting 
from T ~ Tc- This result implies that also in the case of glass-forming liquid the analogously 
defined P quen (q,T) will be trivial when quenching from the paramagnetic (liquid) phase, even 
when quenching down to a T value where RSB occurs (if it does). This is in agreement with 
the behavior reported in |23| for a Lennard- Jones binary mixture starting from T ~ Tmct- 

In glass-forming liquids the overlap among IS is easier to define than the equilibrium 
overlap: its use releases one from the careful consideration of possible crystalline or quasi- 
crystalline configurations. Nevertheless, the quantity P quen (q,T) does not provide any evi- 
dence for the existence or absence of RSB when avoiding the hard task of thermalizing super- 
cooled liquid down to low temperatures. One should also note that in 1RSB models both 
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Figure 1: The equilibrium P(q) (on the left) and the quenched one (on the right) at T = 1.05, 
a value of the temperature slightly higher than Tc — 1, for different values of the system size. 



P(q,T) and P quen (q,T) are expected to be trivial also in the region Tmct > T > Tk (apart 
from finite size effects), just because of the very large number ~ e NYl of 'valleys' (and corre- 
sponding IS) which are almost all orthogonal with zero overlap in the thermodynamic 



limit. 

After clarifying this point, we note that a more careful analysis of the data shown in 
[Fig. |] (right)] suggests the presence of a weak RSB, as already observed in P quen (q,T = oo) 
|^| . In the whole paramagnetic phase, the 'quenched' probability distribution of the overlap 
approaches its delta function limit 5(q) for iV — > oo much slower than the equilibrium one. 
To quantify the iV-dependence, we introduce the function 

/(?) = - lim ^ In [P N , q uen{qj\ ■ (5) 

The replica symmetry is broken in a weak sense if f(q) is zero in an extended region / 
though P qU en(q) is a 5-function in the thermodynamic limit. This implies that PN,quen(q) is 
going to zero slower than exponentially in this region and that therefore, by adding to the 
Hamiltonian a quantity of order iV (for instance by using appropriate boundary conditions), 
one could obtain any given value of the 'quenched' overlap q £ /. Our best numerical evidence 
for such a behavior comes from the study of 

r N (q, T) = — In (^J d qP N , q uen{q, T)^j . (6) 
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Figure 2: r^(q) for the different system sizes at the highest considered temperature T = 3 
(left), at T = 1.8 (center) and at T = 1.2 close to Tc (right). The data plotted by using 
our best numerical estimate v = 0.68 ± 0.05 display a small dependence on N in all the 
paramagnetic phase. 



We get a practically T-independent estimate of the exponent for T > Tc, v = 0.68±0.05. This 
value is compatible with i/ ~ 2/3 obtained from T = oo data in j2"5| . Though one observes 
some deviations at large q, the weak iV-dependence of shown in [Fig. Q and the value of 
v that is significantly smaller than 1, strongly suggest that f(q) is zero on a finite interval /, 
where possibly / = [0, 1], i.e. P que n(q) displays a weak breaking of replica symmetry. 

On the other hand, at temperatures lower than Tc, P quen {q,T) clearly shows the charac- 
teristic behavior corresponding to a full replica symmetry breaking. This is what one would 
expect since RSB was already evident in the P(q) of the configurations we were starting from. 
The qualitative similarities between P(q) and P qU en(q) are remarkable (see [Fig. |||). They 
are present for each disordered configurations: for instance the number of peaks found in a 
one-sample equilibrium Pj(q,T) at a given T is preserved in the corresponding Pj, qU en(q,T) 
too. FRSB features are even more evident when looking at P que n(q)- in particular the pres- 
ence of a continuous plateau between the two self-overlap peaks is clearer in [Fig. |] on the 
right than in [Fig. Q| on the left, since the 'quenched' self overlap takes larger values, though 
it goes to one only for T — > 0. 

Finite size effects play a similar role in the equilibrium and quenched cases: for instance, 
it is clear from [Fig. [| that also the 'quenched' zero overlap probability P quen (0,T) does not 



depend on N in the glassy phase, this being a well known test when looking for FRSB [^3 
To further investigate this point, we plot in [Fig. |[] the ratio of cumulants, 

5(10 = ^(3-^1, (7) 
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(the lowest T-value we have considered) for different system sizes. 
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Figure 4: The Binder parameter from the equilibrium P(q) (left) and from the quenched one 
(right) as a function of temperature for different system sizes. 
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Figure 5: The intersection temperature of the quenched B quen (T) for the different pairs Ni, 
N 2 of considered sizes as a function of Ni + N 2 . 



as a function of T. We recall that this is the usual quantity one calculates in order to locate the 
critical temperature, since finite size scaling predicts that curves for different sizes intersect 
at Tq'- this is the behavior observed in [Fig. |] (left)] (apart from corrections to scaling when 
considering the smaller system sizes). Surprisingly enough, we find that the same kind of finite 
size analysis can be performed on cumulants obtained from probability distribution among 
IS [Fig. H (right)]. Though corrections to scaling are slightly more important in this case, we 
see from [Fig. [|, where the intersection points Tc(Ni, N 2 ) for the different pairs Ni, N 2 of 
considered system sizes as a function of Ni + N 2 are plotted, that one gets the correct T c = 1 
for N h N 2 -> oo. 

We conclude that P que n(q) is an interesting quantity to study. The results we have shown 
could be particularly relevant when looking at glass-forming liquids, but the observed behavior 
suggests that this quantity could also help in further clarifying the glassy phase properties of 
finite dimensional realistic spin glasses (this is a very long-standing subject, see for instance 

mm)- 
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3.2 Correlations between equilibrium configurations and IS 



After noting the similarities between P(q) and P quen (q) one expects the presence of strong 
correlations between equilibrium configurations and corresponding IS, particularly in the low- 
temperature phase. To quantify them, we measure the probability distribution of the overlaps 
q q t between each energy minimum and the configuration from which it is obtained, P q t(q,T). 



it- 




Figure 6: The probability distribution of the overlap between equilibrium configuration and 
corresponding IS P qt (q,T) for different sizes in the high temperature region at T = 3 (left), 
at T = Tc = 1 (center) and in the glassy phase at T = 0.65 (right). 



As it is shown in [Fig. we find a Gaussian shaped distribution which goes toward a 
5-function in the thermodynamic limit both in the paramagnetic and in the glassy phase. We 
observe that below T c , at variance with P(q) and P quen (q), this probability distribution is a 
self-averaging quantity, which is easily understandable since we are substantially looking at 
overlaps between configurations related to the same state. 

It is intriguing to note that there is no clear evidence for the underlying phase transition 
when looking at this quantity. The mean value q q t(T) (see [Fig. 0) is increasing when decreas- 
ing the temperature, and \im T ^ q qt (T) = 1. The iV-dependence appears more pronounced in 
the high temperature region, where data are well fitted by the power law 

C 

qN,qt - qoo,qt+ Jj~ ■ (8) 

We show in [Fig. || our estimates for qoo, q t(T) (left), which gives in particular qoo,qt(Tc) ~ 0.4, 
and our best fit using data for the overlap between random initial conditions and corresponding 
minima (right), i.e. lim^oo qqt{T). Also in this case, we get a non-zero value (~ 0.1) in the 
thermodynamic limit. Nevertheless, it should be stressed that we are confined to a relatively 
small range of system sizes, which makes a reliable estimation of the error of goo,?* very hard. 
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Figure 7: q q t(T) for different sizes as a function of temperature. 
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Figure 8: qoo,qt(T) as obtained by fitting data for different system sizes (left) and our best 
fit by using data on the overlap between random initial conditions and corresponding found 
minima (right), which give limr^oo qoo, q t(T) ~ 0.1. 
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It would be interesting to understand how much these correlations vary when looking at 
different models. In agreement with our results, in | 26fl , for the considered volumes it was 



quoted q qt ~ 0.4 for the SK slightly above T c , to be compared with the higher value q qt ~ 0.94 
found for the 1RSB model (ROM) at a temperature higher than the Mode Coupling T M ct, 
which suggests the presence of stronger correlations between equilibrium configurations and 
corresponding IS in the glass-forming liquid case, and therefore a P quen (q,T) with a behavior 
even closer to the one at the equilibrium. 

3.3 The IS energy 

In [Fig. |9| we show the equilibrium energy e(T), as a function of the temperature T (only 
weakly depending on the system size). In [Fig. |K] (left)] we present data on the IS energy 
£quen{T), i.e. the mean energy of the minima accessible from equilibrium configurations at 
a given temperature and weighted with the Boltzmann factor of the corresponding basin. 
The behavior of this quantity changes abruptly from a nearly T-independent value (in the 
high T regime) to the approximately oc T decreasing of the low temperature region, where 
the IS energy continuously goes towards the ground-state value (which is analytically known, 
e = —0.7633). Correspondingly the derivative de quen (T) / dT , which is plotted in [Fig. [10] 
(right)], displays a maximum and takes very small high-T values (note the logarithmic scale). 

Data on the position T max (N) of the maximum of de quen (T)/dT as a function of iV are 
shown in the insert and give evidence for limjy^oo T max (N) = 1 = Tc- Our finite size analysis 
therefore confirms [25| that in the iV — > oo limit e quen (T) takes the constant threshold value 



2^] eth for T > Tc- The behavior in this region agrees well with the power law 

e N = e th + C/N a , (9) 

giving a constant e#, = — 0.7145±0.004, which is our best numerical estimates for the threshold 
energy, down to T ~ 1.1 near to the critical temperature (and still compatible with this value 
also at Tc). We note that this estimates is in perfect agreement with the e t u — —0.715 quoted 
|5) that was obtained by fitting data on IS reached from random initial conditions by a 



in 



sequential quenching procedure (this value could depend on the considered zero-temperature 
dynamics). The exponent a is slightly increasing when going to lower temperatures and varies 
between a = 0.34 ± 0.04 at T = 3 and a = 0.43 ± 0.06 at T = T C - 

It is interesting to stress that finite size corrections to the asymptotic behavior look very 
important, as shown by the small a value we have found. Correspondingly for all the consid- 
ered sizes (up to the quite large volume N = 1024) the IS energy becomes roughly constant 
only at temperatures definitely higher than Tc- A similar behavior with strong finite size 
corrections is found in |26[ for the considered 1RSB model at T > Tmct too. 
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Figure 9: The equilibrium energy as a function of T for the different system sizes. 
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Figure 10: The quenched energy (left) and its derivative respect to the temperature (right) 
as a function of T for the different system sizes N. In the insert on the right, the temperature 
corresponding to the maximum of de quen /dT is plotted as a function of 1/N. 
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Conclusions 



We have presented numerical results about the properties of energy minima in the SK model. 
The probability distribution of the overlap between IS weighted with the Boltzmann factor of 
the corresponding basin at temperature T, P quen (q,T), turns out to be qualitatively similar 
to the equilibrium overlap distribution P(q, T) at the same temperature T. We found a 
trivially-shaped P quen (q,T) in the whole paramagnetic (T > T c = 1) phase, whereas the 
FRSB behavior characteristic of the glassy phase is evident from data on P quen (q,T) only 
when looking at energy minima obtained from equilibrium configurations at temperature 
definitely lower than Tc- A finite size analysis of the Binder parameter for P quen (q,T) gives 
the same estimate of the critical temperature as the usual (equilibrium) one. 

These results can be particularly relevant for glass-forming liquids, where the overlap is 
more precisely definable between IS [[23|, but they also imply that numerical evidence for 
replica symmetry breaking can not be extracted from data on P quen (q,T), which have been 
obtained from configurations at equilibrium at temperatures above the possible glass transition 
temperature (in that case T# < Tmct)- 

The analysis of finite size corrections to the IS energy confirms |26) that it is approaching 
the expected thermodynamic limit behavior ejs = const = e t h for T > T c , whereas it appears 
continuously decreasing for T <T C . 
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